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Introduction: 

The  CT  concept  (product  of  disinfectant  concentration  and  characteristic 
contact  time)  is  currently  used  to  demonstrate  compliance  with  disinfection  requirements  for 
Giardia  lamblia  ( G .  lamblia )  and  viruses  under  the  Surface  Water  Treatment  Rule  (SWTR).1 
Minimum  CT  requirements  include  relatively  large  safety  factors  to  account  for  possible 
deviations  from  actual  disinfection  efficiencies  achieved  in  full-scale  contactors.  The  application 
of  this  conservative  regulatory  approach  for  Cryptosporidium  parvum  (C.  parvum )  might  result 
in  unrealistic  disinfection  requirements  under  the  Enhanced  SWTR2  due  to  the  much  stronger 
resistance  of  this  protozoan  parasite  to  inactivation  by  all  chemical  disinfectants  used  in  drinking 
water  applications.  There  is  a  need  for  the  development  of  approaches  that  could  provide  a  more 
accurate  assessment  of  actual  inactivation  efficiency  achieved  in  disinfection  contactors. 

The  main  objective  of  this  study  is  to  develop  and  apply  a  mathematical  model  for 
predicting  the  inactivation  of  Cryptosporidium  spp.  (C.  parvum  and  C.  muris )  oocysts  in  ozone 
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bubble-diffuser  contactors.  The  model  is  calibrated  with  semi-batch  kinetic  data,  verified  with 


pilot-scale  inactivation  experiments,  and  used  for  predicting  and  optimizing  full-scale 
disinfection  efficiency. 


Background: 

Waterborne  pathogens  pose  a  substantial  public  health  risk  if  present  in  drinking  water. 
One  common  goal  of  all  potable  water  treatment  plants,  whose  source  is  surface  water,  is  the 
removal  or  inactivation  of  these  pathogens.  Cryptosporidium  is  perhaps  the  most  resistant  to 
disinfection  efforts,  and  therefore  it  has  been  the  subject  of  many  studies  and  well  founded 
concern.  Cryptosporidium  is  classified  as  a  coccidian  protozoan  parasite  of  both  humans  and 
animals.3,4  Initially,  as  many  as  twenty  species  were  named  for  the  host  in  which  they  were 
found;  however,  subsequent  studies  indicated  a  lack  of  host  specificity  thereby  invalidating  many 
of  the  names.  Two  species  affecting  mammals,  Cryptosporidium  parvum  (3.5pm  diameter 
oocyst)  and  Cryptosporidium  muris  (5.8  pm  diameter  oocyst)  have  been  identified.  In  addition, 
two  other  species,  C.  baileyi  and  C.  meleagridis,  are  associated  with  birds.  As  the  name 
suggests,  Cryptosporidium  denotes  an  organism  with  sporozoites  concealed  within  an  oocyst.  As 
it  has  four  aflagellar,  but  motile  sporozoites  with  apical  complexes,  the  genus  Cryptosporidium 
has  been  assigned  to  the  phyllum  Apicomplexa,  class  Sporozoa 

Cryptosporidium  was  first  identified  in  1907  by  Tyzzer  after  being  observed  in  the  gastric 
glands  of  asymptomatic  laboratory  mice.5  It  was  not  until  1976  that  it  was  recognized  as  a 
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human  pathogen  by  Nime.6  Infection  occurs  after  ingestion  and  also  possibly  by  inhalation  of 
the  oocysts.3,4  The  oocysts  are  shed  in  the  feces  of  infected  humans  or  animals.  Transmission 
occurs  by  person  to  person  or  animal  to  person  contact;  ingestion  of  contaminated  food  or  water; 
or  contact  with  contaminated  objects.  Once  inside  the  intestinal  tract,  the  oocysts  reach  the  upper 
small  bowel  where  proteolytic  enzymes,  bile  salts,  and  temperature  enhance  the  excystation  of 
the  four  sporozoites.  The  sporozoites  enter  the  brush  border  surface  epithelium  and  develop  into 
merozoites  capable  of  replicating  either  asexually  or  sexually  beneath  the  cell  membrane  in  the 
brush  border  epithelial  cell  surface.  Sexual  stages  combine  to  form  new  oocysts.  The  new 
oocysts  have  varyings  shell  thicknesses.  Some  of  the  oocysts  generated  with  thinner  than 
average  shells  may  sporulate  and  continue  to  infect  the  same  host,  while  the  others  are  excreted. 
The  cycle  then  continues,  infecting  other  hosts. 

Infection  by  Cryptosporidium  is  called  Cryptosporidiosis.3,4  The  disease  is  characterized 
by  copious,  watery  diarrhea;  abdominal  cramps;  nausea;  vomiting;  anorexia;  weight  loss; 
flatulence;  and  fever.  The  incubation  period  is  2  to  14  days,  with  an  average  of  7  days.  Recent 
studies  suggest  that  the  ID50,  the  dose  required  to  infect  50%  of  the  subjects  tested,  is  132 
oocysts.  The  real  concern  is  that  immunocompromised  individuals  are  at  increased  risk. 
Immunocompromised  individuals  are  young  children,  pregnant  women,  and  people  with  a 
weakened  immune  system  suffering  from  chronic  illnesses  such  as  HIV,  hepatitis,  renal  failure 
and  cancer.  There  is  currently  no  known  effective  therapy  for  cryptosporidiosis.  While  the 
disease  will  typically  run  it’s  course  in  healthy  victims,  cryptosporidiosis  in  immocompromised 
individuals  is  chronic,  progressive,  and  sometimes  fatal.  In  pregnant  women,  cryptosporidiosis 
may  result  in  birth  defects,  miscarriage  or  premature  birth. 


3 


Cryptosporidiosis  is  most  prevalent  in  underdeveloped  countries  where  the  sanitary 
systems  are  lacking  and  drinking  water  is  not  properly  treated.3,4  Numerous  cases  have  also 
occurred  in  the  United  States.  The  most  significant  outbreak  occurred  in  Milwaukee  in  April  of 
1993  and  affected  an  estimated  403,000  people.  Several  immunocompromised  patients  died  in 
what  is  now  considered  the  largest  waterborne  outbreak  in  U.S.  history.  Undoubtedly,  there  have 
been  many  other  outbreaks;  however,  due  to  the  onset  of  symptoms  which  are  very  similar  to 
those  associated  with  a  multitude  of  other  diseases,  many  cases  have  most  likely  either  gone 
unreported,  or  could  not  be  confirmed. 

The  effective  prevention  of  cryptosporidiosis  is  largely  dependent  upon  the  ability  of 
individual  water  treatment  plants  to  either  render  the  oocysts  nonviable,  unable  to  sporulate  and 
reproduce,  through  various  treatment  techniques,  or  to  physically  remove  the  oocysts  from  the 
water  prior  to  distribution.  While  there  has  been  an  abundance  of  recent  research  regarding  the 
chemical  inactivation  of  Cryptosporidium ,  research  is  ongoing.  The  results  currently  appear  to 
indicate  that  the  only  effective  chemical  agents  are  ozone  and  chlorine  dioxide.  Monochloramine 
and  free  chlorine  species  are  generally  considered  ineffective  under  typical  drinking  water 
treatment  conditions7,8,9.  Physical  separation  of  the  oocysts  from  the  water  supply  is  possible 
through  membrane  filtration,  and  to  a  somewhat  lesser  extent  through  conventional  sand 
filtration.  Research  to  quantify  the  effectiveness  of  membrane  filtration,  as  well  as  to  identify  a 
surrogate  indicator  that  would  be  used  in  lieu  of  oocysts,  is  currently  in  progress.  Numerous 
studies  have  documented  the  resistance  of  oocysts  to  extreme  environmental  conditions  for 
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prolonged  periods  of  time  adding  to  the  complexity  and  task  of  achieving  inactivation  to  a  degree 
suitable  to  prevent  human  infection.10,11’12 

A  substantial  problem  with  modem  treatment  techniques  is  the  inability  for  treatment 
plant  operators  to  reliably  measure  viable  oocysts  in  the  drinking  water  they  produce.  There  are 
a  few  techniques,  flow  cytometry  and  epifluorescent  microscopy  for  example,  available  for 
counting  oocysts  in  a  given  volume  of  water;  however,  there  is  no  current  technique  for  a  plant 
operator  to  count  oocysts  leaving  the  plant  and  entering  the  water  distribution  system  and 
distinguish  between  viable  and  inactivated  oocysts.13  Experimentally,  there  are  currently  two 
methods  of  assessing  inactivation  of  viable  Cryptosporidium  spp.  oocysts  in  water  treatment 
systems. 

The  first  viability  assessment  method  is  in-vitro  excystation  which  involves  exposing  the 
oocysts  to  conditions  that  simulate  the  gastrointestinal  environment  of  the  host.  Under  the  proper 
conditions,  presumably  only  the  oocysts  which  are  still  viable  will  excyst  in  accordance  with 
their  life  cycle  and  release  their  four  sporozoites.  The  technique  is  decribed  in  detail  by 
Woodmansee  (1987). 14  Until  a  recent  contribution  by  Rennecker  et  al.  (1997)15,  assessment  of 
inactivation  by  in-vitro  excystation  was  not  very  accurate  due  to  the  inability  to  differentiate 
between  oocysts  that  were  viable  and  had  excysted,  and  those  that  appeared  excysted  but  had 
really  suffered  damage  due  to  over-exposure  to  ozone.  After  excystation  is  promoted,  the 
viability  efficiency  can  be  accurately  and  reproducibly  determined  by  counting  the  excysted 
sporozoites  instead  of  the  empty  shells. 
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The  second  method  is  in-vivo  or  animal  infectivity  which  involves  the  inoculation  of 
animals,  usually  mice,  with  oocyst  suspensions.  The  ID50  is  determined  using  the  Spearman- 
Karber  method  described  by  Finny  (1978). 16  Inactivation  is  then  determined  by  calculating  the 
difference  between  the  infection  dose  for  exposed  oocysts  and  an  unexposed  control.  It  is 
imperative  that  cross  contamination  of  one  animal  by  another  does  not  occur,  otherwise  the  data 
will  be  misinterpreted.  This  technique  can  be  very  time  consuming  and  expensive,  and  can 
additionally  be  the  subject  of  protest  by  animal  activist  groups. 

Both  methods,  however,  can  be  subject  to  error  from  many  sources.  Studies  have 
indicated  that  different  sample  oocyst  cleaning  and  preparation  procedures  may  influence  their 
resistance  and  viability  thereby  altering  data  sets  run  under  the  same  conditions.17  It  is  therefore 
potentially  very  important  that  experiments  are  conducted  with  the  same  batch  of  oocysts. 

In  accordance  with  the  Safe  Drinking  Water  Act  (SDWA),  the  U.S.  Environmental 
Protection  Agency  (U.S.  EPA)  has  proposed  the  Enhanced  Surface  Water  Treatment  Rule 
(ESWTR)  (U.S.  EPA,  1994)  which  will  require  that  water  utilities  provide  average  treatment 
efficiencies  ranging  from  99  to  99.9999%  (2  to  6  logs)  for  C.  parvum  oocysts  through  filtration 
and  disinfection.  The  impact  to  water  utilities  will  be  significant  in  that  they  will  be  forced  to 
determine  what  modifications  and  additions  to  their  water  treatment  plant  are  necessary,  and  then 
further  demonstrate  compliance.  Of  additional  significance  is  that  the  U.S.  EPA  will  have  to 
provide  technical  guidance  and  develop  a  method  of  determining  compliance  with  the  proposed 
standard.  To  compound  matters,  the  simplest  and  more  accurate  method  is  only  realistically 
capable  of  assessing  inactivation  efficiencies  up  to  99.9%  (3  logs).  While  the  human  benefit  is 
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obvious,  there  does  not  presently  exist  a  method  of  measuring  the  exact  inactivation  achieved  by 
disinfection,  so  the  implementation  and  enforcement  of  the  proposed  standard  are  inhibited  by 
technological  limitations  at  this  point  in  time.  The  U.S.  EPA  unfortunately  appears  caught 
between  the  desire  to  protect  the  health  of  the  population  and  the  practicality  of  applying  and 
enforcing  needed  regulations. 

The  drinking  water  industry  is  thus  plagued  with  a  human  pathogen,  Cryptosporidium 
spp.,  which  appears  to  be  present  in  most  natural  waters,  happens  to  be  very  resistant  to  modem 
treatment  techniques  and  is  extremely  difficult  to  measure.  The  federal  agency  charged  with, 
among  many  other  tasks,  the  protection  of  the  health  and  well  being  of  the  population,  has 
encountered  a  very  delicate  situation.  They  must  promulgate  standards  which  will  ensure  the 
safety  of  consumers,  yet  ensure  that  the  capability  exists  to  implement  and  enforce  the  standards. 
It  is  at  this  juncture  that  many  new  and  exciting  research  efforts  have  originated  in  hopes  of 
expanding  the  envelop  of  our  understanding  of  Cryptosporidium  and  it’s  effective  inactivation. 

A  model  accurately  predicting  oocyst  inactivation  could  become  a  great  tool  for  both  regulators 
as  well  as  plant  designing  engineers  and  operators. 


Mathematical  Model  Development: 

The  use  of  ozone  as  a  chemical  disinfectant  is  becoming  more  common  the  drinking 
water  industry.  Recent  research  efforts  have  indicated  that  ozone  is  currently  the  most  effective 
and  safe  chemical  disinfectant  for  pathogens,  requiring  substantially  less  contact  time  to  achieve 


inactivation  than  chlorine,  monochloramine  and  chlorine  dioxide.  The  most  common  method  of 
applying  ozone  is  by  diffusing  it  into  water  inside  bubble-diffuser  contactors  of  varying 
configurations.  Bench  scale  experiments  performed  with  batch  and  semi-batch  reactors  are  being 
used  to  determine  contact  time  and  concentration  required  in  order  to  achieve  a  specific  level  of 
inactivation  for  a  pathogen  of  interest.  There  is  currently  no  tool  available,  however,  to  measure 
the  exact  exposure  time  accrued  in  flow-through  bubble-diffuser  ozone  contactors. 

Currently,  the  CT  concept  is  used  to  demonstrate  compliance  with  disinfection 
requirements  under  the  SWTR.  Typically,  a  t10  hydraulic  residence  time  is  used  to  predict  the 
contact  time.  This  contact  time  will  be  conservative  for  most  reactor  designs  except  those 
approaching  perfectly  mixed  (CSTR)  conditions.  The  t10  hydraulic  residence  time  is  calculated 
from  a  non-reactive  tracer  test  conducted  on  the  contactor.  It  is  defined  as  the  period  of  time  that 
it  takes  for  10%  of  the  mass  of  a  pulse  input  to  reach  the  effluent,  or  alternatively  for  the  effluent 
concentration  to  increase  to  10%  of  the  net  concentration  of  a  step  input.  The  tracer  test  provides 
an  indication  of  the  amount  of  time  that  90%  of  the  oocysts  entering  the  contactor  at  a  given  time 
will  physically  remain  in  the  contactor.  The  t10  time  is  based  on  the  mixing  and  dispersive 
characteristics  of  the  ozone  diffuser  and  chamber  configuration.  The  t,0  time  is  then  multiplied 
by  the  average  dissolved  ozone  concentration  in  a  given  chamber  resulting  in  total  CT.  If 
research  indicates  that  5  minutes  of  “contact  time”  are  required  per  mg/1  of  a  known  pathogen  to 
achieve  3  logs  (99.9%)  of  inactivation,  the  plant  operator  must  ensure  that  the  t10  time,  which  is 
fixed  for  a  given  contactor,  multiplied  by  the  average  ozone  concentration,  results  in  5  min-mg/1. 
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The  plant  operator  is  then  left  with  the  option  of  altering  the  quantity  of  ozone  injected 
into  the  chamber  as  the  sole  means  of  achieving  a  desired  level  of  inactivation.  This  approach  is 
conservative  and  may  additionally  result  in  unrealistic  disinfection  requirements  specifically  for 
C.  parvum  under  the  ESWTR  due  to  the  much  stronger  resistance  of  this  pathogen  to  inactivation 
by  all  chemical  disinfectants  used  in  drinking  water  applications.  Equally  challenging  is  the 
design  phase  of  an  ozone  contactor.  The  designer  can  modify  height,  diameter,  diffuser  type,  add 
additional  chambers,  split  the  ozone  input  over  two  or  more  different  chambers,  specify  counter 
or  co-current  flow,  and  alter  the  liquid  and  gas  flow  rates.  While  the  designer  has  software 
available  to  predict  impacts  on  cost  and  physical  performance,  there  is  no  existing  tool  with 
which  performance  related  to  inactivation  of  pathogens  can  be  evaluated. 

The  intent  of  the  model  developed  in  this  study  is  to  provide  a  basis  for  making  sound 
decisions  about  the  design  and  operation  of  ozone  contactors  in  an  effort  to  achieve  a  desired 
goal  of  pathogen  inactivation.  The  model  is  formulated  from  mass  balance  principles 
incorporating  contactor  hydrodynamic  information  obtained  from  tracer  tests,  experimentally 
measured  pathogen  inactivation,  ozone  decomposition  and  mass  transfer  kinetics.  The 
hydrodynamics  inside  the  contactor  chambers  are  assumed  to  approach  ideal  axial  dispersion 
reactor  (ADR)  conditions.  The  mass  balance  equations  provide  an  accurate  prediction  of  the 
dissolved  ozone  concentration  profile  throughout  a  given  chamber.  The  designer  may  select 
either  counter  or  co-current  flow  on  which  to  base  the  contactor.  Counter  current  as  implied,  is  a 
chamber  in  which  ozone  gas  enters  from  the  bottom  and  the  water  enters  from  the  top.  Both 
ozone  gas  and  liquid  influent  enter  from  the  bottom  in  a  co-current  chamber.  Any  number  of 
additional  chambers  can  represented.  Ozone  may  be  applied  in  any  chamber.  If  no  ozone  is 
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applied  in  a  given  chamber,  the  chamber  simply  becomes  reactive,  in  which  the  dissolved  ozone 
exiting  the  previous  chamber  is  allowed  to  continue  decomposing,  meanwhile  providing 
additional  contact  time  with  the  microorganisms  with  perhaps  less  back  mixing. 

The  equations  representing  the  calculations  of  all  basic  parameters  required  for  any 
configuration  are  given  in  Appendix  I.  Using  the  model,  an  operator  may  input  the  physical 
dimensions  and  flow  rates  of  his  or  her  current  plant  including  chamber  height  and  diameter, 
water  temperature  and  density,  liquid  and  gas  flow  rates,  concentration  of  ozone  applied,  and  the 
dispersion  characteristics  (dispersion  number)  of  the  contactor.  The  dispersion  number  is  the 
ratio  of  liquid  diffusivity  to  the  liquid  velocity  and  chamber  height  and  represents  the  extent  of 
physical  mixing  within  the  contactor.  If  the  dispersion  number  is  unknown,  it  may  be  estimated 
from  tracer  data  using  an  iterative  process  employing  the  Thomas  Method.18  This  technique  is 
demonstrated  later  when  a  contactor  owned  by  the  U.S.  EPA  Research  Laboratory  in  Cincinnati 
is  modeled  in  order  to  compare  the  results  and  accuracy  of  this  model’s  predictions. 

The  designer  is  at  liberty  to  create  virtually  any  dimensions  and  parameters  in  order  to 
observe  the  impact  on  inactivation  achieved.  Both  users,  the  designer  as  well  as  the  operator, 
must  be  careful  to  account  for  the  ozone  demand  in  the  influent  water,  a  parameter  which  may 
easily  vary  with  season,  rain  events,  or  organic  matter  concentration.  Every  body  of  water  will 
have  a  natural  ozone  demand  based  primarily  on  the  organic  content  of  the  water.  Designers 
must  be  sure  to  use  an  ozone  decomposition  rate  applicable  to  the  body  of  water  under  worse 
conditions  under  which  the  contactor  is  to  operate.  Last,  any  pathogen  can  be  modeled,  provided 
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its  inactivation  kinetics  are  known.  The  user  is  referred  to  previously  published  inactivation  rates 
for  a  specific  pathogen  of  interest. 

The  equations  resulting  from  a  mass  balance  on  a  counter-current  chamber  configuration 
are  derived  in  detail  in  Appendix  II.  The  overall  mass  balance  equation  on  the  liquid  side  is 
represented  by: 


°C±  =  E  <?2Q 

a  1  ax2 


s  cK 


+  KL-a-(cl-CL)-KD-CL 


(1) 


The  equation  accounts  for  turbulent  diffusivity,  advection,  mass  transfer  from  the  gas  to 
the  liquid  phase  governed  by  the  Two  Film  Mechanism  and  Henry’s  Law,  and  the  decomposition 
of  the  ozone  in  the  liquid  phase.  The  liquid  volume  fraction,  e,  can  be  taken  as  approaching 
unity  for  typical  drinking  water  conditions,  and  the  corresponding  steady-state  equation  can  be 
partially  non-dimensionalized  to: 


dC L_ 
dZ 


-Nd-Cl=  0 


(2) 


The  overall  mass  balance  equation  on  the  gas  side  is  represented  by: 


_F  Up  0CG 

a  G  ax2  i -e  ax 
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r 

d  W; 


(3) 
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The  equation  accounts  for  turbulent  diffusivity,  advection,  mass  transfer  from  the  gas  to 
the  liquid  phase,  and  the  decomposition  of  the  ozone  in  the  gas  phase.  Making  the  assumptions 
that  diffusivity  and  reactivity  of  ozone  in  the  gaseous  phase  are  negligible,  assuming  steady  state, 
taking  the  liquid  volume  fraction,  s,  as  approaching  unity,  and  non-dimensionalizing  each  term, 
the  equation  can  be  reduced  to: 


dZ 


S 


(Cr. 


-C, 


m 


=  0 


(4) 


Equations  2  and  4  can  be  solved  after  applying  Danckwert’s19  boundary  conditions  and 
the  known  ozone  gas  input  at  the  contactor  bottom.  The  solution  provides  the  following 
expression  for  the  dissolved  ozone  profile  through  the  contactor: 


(  „  \  f  r.  \  f  c 

CL(z)=  \EXP{Axz)-—  AxEXP(Axz)  ax  +  EXP  (/l2Z)-  -^—h2EXP  (d2Z)  a2  +  EXP  (A3Z)-  -—A3EXP  (/l3z) 
V  NL  )  V  NL  J  \ 

(5) 


This  profile  will  be  important  and  used  later  in  predicting  pathogen  inactivation 
throughout  a  contactor  chamber.  Appendix  III  indicates  the  development  of  equations  resulting 
from  a  mass  balance  on  a  co-current  chamber  configuration.  Similarly,  Appendix  IV  indicates 
the  development  of  equations  resulting  from  a  mass  balance  on  a  reactive  chamber.  The  equation 
representing  a  reactive  chamber  is  simplified  since  there  is  no  mass  balance  on  the  gas.  The 
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mass  balance  on  the  liquid  will  contain  only  advection  and  reaction  terms  if  plug  flow  conditions 
are  assumed.  At  this  juncture,  it  should  be  noted  that  the  reactive  chamber  mass  balance 
equations  can  additionally  be  utilized  in  the  predicted  performance  of  a  contacting  basin,  such  as 
a  chlorine  contactor,  where  the  introduction  of  a  gas  is  not  an  issue. 

With  a  complete  profile  of  dissolved  ozone  throughout  every  chamber  in  an  ozone 
contactor,  a  mass  balance  on  the  pathogen  of  interest  can  be  performed.  The  resulting  equations 
are  derived  in  detail  and  presented  in  Appendix  V.  The  overall  mass  balance  on  the  pathogen  is 
represented  by: 


=  e  <rN 
dt~  L'  8X 


■  —  -Kn-N-Cl 

ax. 


(6) 


The  equation  accounts  for  turbulent  diffusivity,  advection,  and  the  inactivation  reaction  of 
the  pathogen  with  ozone.  If  steady  state  is  assumed,  and  the  equation  is  partially  non- 
dimensionalized,  the  equation  can  be  reduced  to: 


d- 


d2N  dN 


dZ1  dZ 


- Nn  •  N -CL  =  0 


(7) 


The  dependence  of  pathogen  inactivation  on  dissolved  ozone  concentration  is  readily 
apparent.  Since  the  concentration  of  dissolved  ozone  is  not  constant  throughout  the  chamber,  an 
iterative  finite  element  approach  is  taken.  The  resulting  tri-diagonal  matrix  can  be  solved  by  the 
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Thomas  Method.  Finally,  an  expression  is  derived  for  the  inactivation  of  the  pathogen  which 
accounts  for  the  history  of  exposure  through  a  changing  ozone  profile. 

The  reaction  term  is  represented  by  the  Chick-Watson  Law  of  pseudo-first  order 
inactivation  kinetics.20  The  Chick-Watson  Law  can  be  applied  in  the  case  of  C.  muris,  as  its’ 
inactivation  conforms  to  that  of  a  pseudo-first  order  reaction.  However,  the  Chick-Watson  Law 
may  not  be  applicable  to  all  pathogens,  in  which  case,  the  correct  mathematical  or  empirical 
relationship  should  be  used.  Chen  (1998)21  devised  an  ingenious  method  to  account  for  the 
inactivation  of  a  pathogen  which  indicates  a  “shoulder”  profile.  An  example  is  C.  parvum  which, 
as  indicated  by  Rennecker  et  al.  (1997)15,  displays  a  resistance  to  ozone  until  a  certain  amount  of 
contact  time  is  exceeded,  then  follows  the  Chick-Watson  first  order  reaction  rate  kinetics.  Chen 
(1998)21  recommends  linear  regression  of  the  inactivation  profile  that  does  conform  with  a  first 
order  reaction  rate  in  order  to  find  a  theoretical  starting  amount  of  pathogens,  N/N0,  higher  than 
100%.  With  the  proper  programming  statements,  the  model  can  represent  any  amount  of 
pathogens  calculated  to  be  greater  than  100%  as  only  100%,  indicating  no  achieved  inactivation. 
The  model  would  then  display  a  “shoulder”,  or  period  of  no  inactivation,  until  the  first  order  rate 
law  would  take  effect  and  begin  indicating  pathogen  amounts  of  less  than  100%.  This  technique 
is  demonstrated  later  when  a  contactor  owned  by  the  U.S.  EPA  Research  Laboratory  in 
Cincinnati  is  modeled  in  order  to  compare  results  and  accuracy  of  this  model’s  predictions  with 
the  experimental  inactivation  of  C.  parvum. 

The  equations  developed  in  Appendices  I  through  V  are  ideal  for  spreadsheet  application 
where  graphs  can  be  used  for  a  visual  representation  of  both  the  dissolved  ozone  and  pathogen 
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inactivation  profiles.  The  spreadsheet  provides  a  user-friendly,  easy  way  of  inputting  and 
altering  design  and  operational  parameters  and  observing  the  inactivation  achieved.  The  model 
can  then  be  used  to  evaluate  the  effectiveness  of  any  configuration  on  the  inactivation  of  a 
pathogen  of  interest.  The  designer  can  optimize  performance  while  the  operator  can  understand 
the  impact  of  altering  operational  parameters.  This  model  combines  the  results  of  extensive 
research  with  the  mathematical  characteristics  of  the  application  and  provides  a  useful  tool  which 
summarizes  anticipated  performance. 


Measuring  Model  Accuracy: 

In  order  to  validate  this  model,  experimental  data  from  Owens  et  al.  (1994)22  was  used. 
Owens  and  co-workers,  employed  by  the  U.S.  Environmental  Protection  Agency  Research 
Laboratory  in  Cincinnati,  Ohio,  have  conducted  several  experiments  with  their  pilot  scale  bubble 
diffuser  depicted  by  Figure  1.  Their  single  chamber  ozone  contactor  measures  2.65  meters  in 
height  and  0.15  meters  in  diameter  and  is  operated  in  the  counter-current  configuration.  It  has 
sampling  ports  at  the  influent  and  effluent,  as  well  as  at  four  intermediate  locations  spaced 
equally  at  0.46  meters  apart.  The  ozone  was  produced  from  oxygen  by  a  model  GL-1  generator 
made  by  PCI,  Inc.  The  concentration  of  ozone  applied  to  the  contactor  and  that  released  in  the 
off-  gas  was  measured  by  ultra  violet  (UV)  light.  The  difference  in  the  two  measurements  gave 
the  amount  of  ozone  transferred  to  the  liquid.  The  transfer  efficiency  was  estimated  to  be  greater 
than  94%.  Additionally,  the  ozone  demand  was  estimated  by  measuring  the  ozone  concentration 
in  the  effluent  at  the  last  sampling  port  and  subtracting  it  from  the  ozone  transferred.  The  liquid 
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Influent 


Figure  1 


Diagram  of  pilot-scale  ozone  contactor 


and  gas  flow  rates  used  were  6.4  L/min  and  0.64  L/min,  respectively.  Step  dose  tracer  studies 
revealed  a  theoretical  mean  residence  time,  Te,  of  7.4  minutes  and  a  T10  time  of  2.27  minutes. 
Filtered  Ohio  River  water  was  used  throughout  all  experiments.  Experimantal  conditions  were  a 
temperature  of  22-25°C,  pH  of  7.5-8.5,  and  a  total  organic  content  (TOC)  of  1.61  to  2.08  mg/L. 


Oocysts  were  obtained  from  various  sources,  Giardia  muris  from  the  University  Hospital 
in  Cleveland  Ohio  from  female  mice,  Cryptosporidium  muris  from  Japan’s  Osaka  University, 
Medical  School  also  from  female  mice,  and  Cryptosporidium  parvum  from  a  locally  infected 
Holstein  calf.  A  suspension  of  1 00  oocysts/mL  was  fed  at  2mL/minute  into  the  filtered  Ohio 
River  water,  just  ahead  of  the  contactor.  Samples  were  then  collected  at  the  various  ports  within 
the  contactor  and  assayed  via  both  in  vitro  excystation  and  animal  infectivity.  They  were  then 
able  to  display  inactivation  profiles  through  the  contactor  from  data  generated  from  the 
excystation  and/or  infectivity  assays. 

Comparing  the  predictions  of  the  mathematical  model  to  the  results  obtained 
experimentally  will  presumedly  validate  the  accuracy  of  the  model.  A  few  parameters  required 
to  run  the  model  were  not  specifically  known,  and  thus  had  to  be  either  derived  or  empirically 
calculated.  First,  while  the  transferred  ozone  concentration  was  determined,  the  actual 
percentage  of  pure  ozone  applied  at  the  bottom  of  the  contactor  was  not  known.  By  using  a 
simple  mass  balance: 

Qg  *  C^gin  "  Cgout)  =  Ql*Q.  ^ 
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and  by  assuming  the  transfer  efficiency  of  “greater  than  94%”  was  97%,  effectively  reducing  the 
maximum  possible  error  in  calculation  to  3%,  CG0UX  could  be  calculated  to  be  0.03  CGIN. 
Knowing  the  transferred  ozone  concentration,  CL,  CGIN  could  be  determined  with  the  following 
relationship  for  each  of  the  three  levels  of  reported  transferred  ozone: 


-GIN  ' 


Ql. 

Qg 


*CL* 


1 

0.97 


(9) 


Additionally,  the  dispersion  number  for  the  contactor  was  unknown.  However,  since  the 
tracer  tests  produced  T0and  T10  times,  the  ratio,  T1O/T0  could  be  calculated  to  be  0.3068.  Using 
the  Thomas  Method,  described  in  detail  in  Appendix  V,  trial  values  for  d,  the  dispersion  number, 
can  be  generated  until  the  resulting  profile  indicates  a  T1O/T0  value  of  0.3068.  Trial  and  error 
results  in  a  dispersion  number  of  0.4424  which  can  be  used  to  model  the  EPA’s  contactor.  The 
Thomas  Method  is  very  useful  and  is  commonly  used  as  a  tool  to  integrate  tracer  concentrations 
over  time  to  arrive  at  a  step  profile. 

Last,  the  model  applies  a  known  ozone  demand  to  predict  inactivation.  Since  ozone 
demand  is  a  temperature  dependent  reaction  term,  the  data  accumulated  by  Owens  et  al.  must  be 
translated  into  a  reaction  coefficient,  KD.  Recall  the  equation  given  by  Appendix  II  for  the  liquid 
side: 
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(10) 


dCj_  =  E  ?CL 
a  L  ax2 


ax 


+  KL-a(CL-CL)-KD-CL 


In  this  case,  since  the  dissolved  ozone  concentration  in  the  contactor  effluent  was  reported,  the 
reaction  coefficient,  KD,  can  be  altered  until  the  proper  resulting  values  of  dissolved  ozone 
effluent  concentration  is  properly  matched  by  the  model.  Table  1  summarizes  all  relevant 
parameters. 

At  this  point,  the  model  has  all  the  parameters  required  to  predict  a  dissolved  ozone 
concentration  profile  throughout  the  contactor.  Additionally,  if  the  kinetic  reaction  term,  KN, 
applicable  to  a  given  pathogen  is  known,  the  model  will  predict  an  inactivation  profile  for  the 
pathogen  throughout  the  contactor.  Inactivation  of  pathogens  are  given  by  the  Chick-Watson 
pseudo-first  order  reaction  rate:  Kn*N*Cl,  where  KN  is  the  second  reaction  rate  constant,  N  is 
the  count  of  active,  viable  oocysts  at  a  given  point  in  the  contactor,  and  CL  is  the  corresponding 
dissolved  ozone  concentration  at  that  point.  C.  muris  conforms  ideally  to  this  pseudo-first  order 
reaction.  Unfortunately,  previous  research  documented  by  Rennecker  et  al.  (1997) 15  only 
accounted  for  the  inactivation  of  C.  muris  through  a  temperature  range  between  5  and  20  °C, 
while  the  experiments  completed  by  Owens  et  al.  (1994)22  were  conducted  at  22  to  25  °C. 
Rennecker  et  al.  (1997)15  suspected  that  above  20  °C,  the  inactivation  coefficient  increased  at  a 
much  greater  rate  than  the  temperature  dependency  indicated  between  5  and  20  °C.  By 
following  the  same  procedures  outlined  in  Rennecker  et  al.  (1997)15,  additional  experiments  were 
conducted  with  C.  muris  which  has  provided  a  temperature  dependent  inactivation  coefficient 
valid  through  30  °C  represented  by  the  following  expression: 
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C  muris  Experimental  Data  (Owens  et  al.  1994) 


Transferred 

Ozone 

(mg/L) 

Contactor 

Effluent 

(mg/L) 

Ozone 

Demand 

(mg/L) 

Reaction 

Constant 

kd 

Log 

Inactivation 

1.03 

0.37 

.416 

0.41 

1.37 

0.32 

1.05 

.867 

0.51 

1.38 

0.34 

1.04 

.789 

0.36 

2.07 

1.22 

.129 

1.20 

2.11 

1.26 

.124 

1.58 

2.76 

1.54 

1.22 

.154 

1.64 

3.31 

1.87 

1.44 

.148 

2.56 

3.72 

2.05 

1.67 

>2.70 

4.63 

2.53 

2.10 

.165 

>2.62 

Temperature:  23.6  +/- 1.6  °C 


C.  parvum  Experimental  Data  (Owens  et  al.  1994) 


Transferred 

Ozone 

(mg/L) 

Contactor 

Effluent 

(mg/L) 

Ozone 

Demand 

(mg/L) 

Reaction 

Constant 

kd 

Log 

Inactivation 

1.58 

0.87 

.262 

1.08 

1.76 

0.88 

0.89 

.198 

0.57 

2.71 

1.46 

1.25 

.163 

1.80 

3.20 

1.69 

1.51 

.173 

2.17 

3.31 

1.81 

1.50 

.157 

1.95 

3.56 

1.71 

1.84 

.2205 

2.67 

3.97 

2.35 

1.46 

.124 

4.17 

2.31 

1.86 

.152 

>2.15 

Temperature:  24.5  +/-  1.6  °C 


Table  1 


Experimental  parameters 


Kn  =  EXP(12.5)  X  EXP 


(  33500  ^ 


,  where  temperature  is  in  °K 


(ID 


8.314xT ) 

When  this  relationship  is  used,  the  model  predicts  inactivation  per  Figure  2.  Data  from  Owens  et 
al.  (1994)22  is  depicted  also  to  illustrate  the  close  match  between  experimental  and  model 
predicted  inactivation. 

Some  pathogens,  C.  parvum  for  example,  may  not  ideally  conform  to  the  Chick-Watson 
first  order  inactivation  rule.  C  parvum  actually  indicates  a  shoulder,  a  period  of  no  inactivation 
with  ozone  exposure  until  a  certain  contact  time  threshold  is  exceeded.  After  reaching  the 
threshold,  the  first  order  inactivation  rate  is  applicable.  Figure  3  indicates  the  inactivation  profile 
for  C.  parvum.  If  the  linear  portion  of  the  inactivation  profile,  that  which  does  comply  with  a 
first  order  rate  law,  is  projected  upwards  to  an  intercept  above  N/No=100,  the  entire  profile  can 
be  modeled  as  first  order,  per  Figure  4.  This  means  that  the  model  can  be  used  to  simulate  a 
starting  pathogen  count  of  higher  than  100%  and  make  use  of  the  Chick-Watson  first  order  law  to 
predict  the  entire  profile.  With  some  logic  built  into  the  programming,  the  model  can  ignore  all 
inactivation,  N/N0,  profile  calculations  predicted  to  be  in  excess  of  100  and  simply  letN/No=100 
for  those  data  points.  In  this  manner,  the  model  can  account  for  a  “shoulder”.  Work  by 
Rennecker  et  al.  (1997)15  has  provided  a  temperature  dependent  inactivation  coefficient  for  C. 
parvum  represented  by: 


KN=  1.2985  X  1015  EXP 


10079.4'' 
T  j 


,  where  temperature  is  in  °K. 


(12) 
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0.1  1  10  100 

C.  muris  Survival  (%) 


Figure  2 


Comparison  of  model  predicted  to  experimentally  determined  inactivation  of  C. 


C.  parvum  Inactivation  "shoulder" 


Contact  Time  (mg  -  min/L) 


•  Data  Points  from  Rennecker  et  al.  (1998) 


Figure  3 

Experimental  data  from  Rennecker  et  al  (1998)  depicting  the  shoulder  trend  observed  with  C. 
parvum.  Shoulder  signifies  a  required  contact  time  prior  to  any  observed  inactivation  beginning 


Extrapolation  of  new  N0  to  Account  for  Shoulder 


Contact  Time  (mg-min/L) 


•  Data  Points  from  Rennecker  et  al.  (1998) 
-  Linear  Regression 


Figure  4 

By  linear  regression,  a  new  N/N0  is  obtained  in  order  to  make  use  of  the  entire  linear  profile 
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When  this  relationship  is  used,  the  model  predicts  inactivation  per  Figure  5.  Data  from  Owens 
et  al.  (1994)22  is  depicted  also  to  illustrate  the  close  match  between  experimental  and  model 
predicted  inactivation. 

As  can  be  observed  from  Figures  2  and  5,  the  model  predicts  inactivation  very  accurately 
as  compared  to  the  inactivation  measured  experimentally  by  Owens  et  al.  (1994)22. 

Discrepancies  can  be  attributed  to  many  factors  inherent  in  any  experiment.  First,  the  exact 
temperature  of  the  Ohio  River  water  was  not  known,  but  rather  a  range  was  provided.  This  is  not 
viewed  as  a  significant  error  because  the  model  can  simulate  the  relatively  small  range  and 
predict  a  corresponding  “range”  of  inactivation  which,  when  done,  would  be  rather  small. 

Second,  the  ozone  demand  for  a  given  body  of  water  will  vary  significantly  from  month  to  month 
and  following  rain  events.  The  experiments  run  by  Owens  et  al.  (1994)22  were  conducted  over  a 
period  of  a  few  months,  necessitated  by  the  scarcity  of  available  oocysts  and  the  physical 
duration  of  excystation  and  infectivity  tests.  Third,  an  empirical  dispersion  number  was 
generated,  in  lieu  of  knowing  the  specific  contactor  dispersive  characteristics.  Fourth,  the  tests 
conducted  by  Owens  et  al.  (1994)22  can  measure  with  a  certain  degree  of  accuracy  up  to  2  logs  of 
inactivation  when  using  the  excystation  method,  however,  beyond  that  point  lose  reliability. 
Effectively,  to  measure  for  example  3  logs  of  inactivation  (99.9%)  by  excystation,  a  count  of 
8000  oocysts  is  required  to  obtain  a  statistically  accurate  answer.  Last,  removing  samples  from 
ports  within  a  contactor  and  ensuring  that  quenching  occurs  immediately  is  critical.  Otherwise, 
residual  dissolved  ozone  may  continue  to  effect  the  oocysts,  even  briefly,  resulting  in  an 


24 


C.  pa/vumSurvival  (%) 


ADR  Model  Fitting  of  C.  parvum  Oocyst  Inactivation  Data 


•  Data  by  Owens  et  al.  (1994) 

°  ADR  Model  Fitting 

—  Linear  Regression 


Figure  5 


Comparison  of  model  predicted  to  experimentally  determined  inactivation  of  C.  parvum 


overestimation  of  inactivation.  Additionally,  if  the  sampling  port  comes  off  a  pipe  even  a  few 
inches  long,  plug  flow  conditions  with  high  dissolved  ozone  concentrations,  especially  at  the 
bottom  sampling  ports,  may  likewise  induce  an  overestimation  of  inactivation. 

Taking  into  account  all  the  possible  inaccuracies,  the  model  predicts  extremely  well 
compared  to  the  experimental  data  of  Owens  et  al.  (1994)22,  which  to  date  is  really  the  only 
known  available  data  with  which  to  verify  the  model.  The  observed  margin  of  difference 
between  the  experimental  and  model  predictions  is  on  the  order  of  magnitude  of  less  than  1  %. 
The  model’s  predictions  are  certainly  within  the  range  of  accuracy  suitable  to  forego  expensive 
and  time  consuming  experimental  tests,  which  are  themselves,  inherently  subject  to  error. 


Model  Applications: 

With  the  model  set  up  on  a  spreadsheet,  the  applications  become  endless.  The  designer 
of  a  bubble-diffuser  ozone  contactor  is  at  liberty  to  alter  almost  any  parameter  and  observe  the 
resultant  effect  on  inactivation  predicted.  The  plant  operator  who  is  managing  an  existing 
contactor  in  operation  is  somewhat  restricted,  however,  they  can  still  significantly  alter  the 
inactivation  by  a  few  very  important  parameters.  The  usefulness  of  the  model  in  optimizing 
plant  efficiency  can  be  observed  by  a  few  examples. 


To  begin  with,  the  bubble-diffuser  can  be  designed  as  a  simple  one  chamber  system, 
however,  it  may  be  configured  for  counter-current  or  co-current.  Using  the  dimensions  and 
operating  parameters  of  the  contactor  previously  described,  the  model  was  run  for  various 
applied  ozone  concentrations.  As  depicted  by  Figure  6,  in  each  case,  higher  inactivation  was 
achieved  for  the  co-current  configuration.  Changes  to  the  height  and  diameter  of  the  chamber 
did  not  change  the  result.  The  verdict  is  that  if  a  single  chamber  model  is  desired,  it  should  be 
configured  for  co-current  operation. 

If  multiple  chambers  are  desired,  the  results  are  interestingly  not  as  predictable.  For 
chambers  all  of  equal  diameter,  the  counter-current  configuration  achieves  higher  inactivation 
compared  to  the  co-current  configuration  at  varying  concentrations  of  applied  ozone  per  Figure 
7.  When  the  chamber  diameters  are  altered  such  that  chambers  one  and  three  are  larger,  15  cm 
for  example,  and  chambers  two,  four,  five  and  six  are  smaller,  5  cm  for  example,  the  counter- 
current  configuration  can  be  optimized  and  achieves  higher  inactivation  for  the  same 
concentrations  of  applied  ozone.  Figure  8  indicates  the  observed  inactivation  corresponding  to 
the  two  configurations.  By  viewing  the  concentration  of  dissolved  ozone  in  the  effluent  of  the 
sixth  chamber,  it  can  be  seen  that  there  is  a  slightly  higher  concentration  of  ozone  in  the  effluent 
of  the  co-current  configuration.  It  would  seem  that  a  lower  inactivation  is  achieved  potentially 
due  to  more  of  the  ozone  not  being  utilized.  The  counter-current  configuration  results  in  a 
slightly  smaller  dissolved  ozone  concentration  in  the  effluent  indicating  more  efficient  use  of  the 
ozone  towards  inactivation. 
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N/N0  (LOG  Inactivation  Achieved) 


One-Chamber  Model  Comparison 


Applied  Ozone  Concentration  (mg/L) 


Counter-Current  Model  Predictions 
o  Co-Current  Model  Predictions 


Figure  6 

Comparison  of  inactivation  effectiveness  of  one-chamber  design  configurations 


Multi-Chamber  Model  Comparison 
(Constant  Chamber  Diameters) 


Applied  Ozone  Concentration  (mg/L) 


« —  Counter-Current  Model  Predictions 
o  Co-Current  Model  Predictions 


Figure  7 


Comparison  of  inactivation  effectiveness 
configurations 


of  multi-chamber,  constant  diameter  design 
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Multi-Chamber  Model  Comparison 
(Varying  Chamber  Diameters) 


Applied  Ozone  Concentration  (mg/L) 


Counter-Current  Model  Predictions 
o  Co-Current  Model  Predictions 


Figure  8 


Comparison  of  inactivation  effectiveness  of  multi-chamber,  variable  diameter  design 
configurations 
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As  is  frequently  the  case,  many  multi-chamber  bubble  diffusers  permit  split  ozone  input. 
To  accomodate  and  predict  the  efficiency  of  these  designs,  the  model  allows  for  the  input  of 
ozone  in  chambers  one  and  three.  The  interpretation  of  the  results  in  the  comparison  of  both 
configurations  is  not  as  obvious.  For  varying  concentrations  of  ozone  applied,  and  varying 
chamber  diameters,  the  co-current  model  achieved  highest  inactivation  consistently  when  100% 
of  the  ozone  was  applied  in  the  first  chamber.  While  not  performing  as  well  as  the  co-current 
configuration,  the  counter-current  configuration  indicated  that  inactivation  can  be  maximized 
with  a  75/25  %  split  ozone  input  between  chambers  one  and  three  if  dispersion  was  assumed  in 
the  reactive  chambers.  When  assuming  no  dispersion  occurs  in  the  reactive  chambers,  plug  flow 
conditions  prevail  and  inactivation  can  be  maximized  with  a  25/75  %  split  ozone  input  between 
chambers  one  and  three.  This  only  held  true  in  the  case  where  the  chamber  diameters  were  equal 
and  there  was  an  indication  that  the  ozone  was  depleted  too  much  when  100%  was  applied  in  the 
first  chamber  to  make  effective  use  of  the  last  chamber.  Again,  the  results  indicated  that  the  level 
of  inactivation  was  related  to  the  efficient  use  of  dissolved  ozone  so  that  neither  high  nor  very 
low  concentrations  were  observed  in  the  effluent  from  the  final  chamber.  A  high  final  chamber 
effluent  concentration  indicated  a  waste  of  ozone  while  a  low  effluent  concentration  indicated  a 
waste  of  a  chamber. 

Seasonal  fluctuations  in  temperature,  as  may  be  imagined,  could  pose  a  substantial 
impact  to  plant  efficiency.  The  exact  relationship  of  temperature  dependence  is  another  factor 
which  is  not  inherently  obvious  and  the  model  may  be  of  great  assistance.  Using  data  and 
parameters  reported  by  Owens  et  al.  (1994)22,  it  may  be  observed  that  an  increase  in  temperature 
results  in  greater  inactivation  of  C.  parvum  per  Figure  9.  This  is  not  a  result  which  should  be 
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taken  as  a  conclusive  trend  in  every  case.  It  is  possible  that  specific  sources  of  drinking  water 
may  have  large  seasonal  variations  in  organic  content.  If  the  ozone  reaction  or  degradation 
coefficient  is  adjusted  to  reflect  increasing  organic  content  with  higher,  summer  temperatures,  it 
is  possible  that  ozone  demand,  if  excessive,  could  actually  result  in  lower  inactivation.  The 
conclusion  is  thus  not  as  predictable  without  the  help  of  the  model  and  accounting  for  the 
conditions  inherent  to  the  specific  water  source.  The  designer  and  operator  must  account  for 
ozone  demand  in  the  influent  water  before  making  decisions  based  on  temperature  alone.  With 
the  model,  the  operator  may  make  an  intelligent  decision  to  either  increase  or  reduce  the  amount 
of  applied  ozone  as  the  seasonal  temperature  and  ozone  demand  changes  in  order  to  consistently 
meet  inactivation  requirements. 

The  results  obtained  when  modifying  a  few  parameters  are  interesting  and  not  readily 
predictable.  They  indicate  that  the  operator,  while  physical  dimensions  and  parameters  can  not 
be  altered,  has  a  wide  range  of  flexibility  and  can  indeed  optimize  the  plant’s  efficiency 
regarding  pathogen  inactivation.  A  designer  has  unlimited  flexibility,  and  if  specifically  aware 
of  the  influent  characteristics,  temperature  and  ozone  demand,  can  make  very  wise  engineering 
decisions  regarding  plant  configuration,  number  of  chambers,  applied  ozone,  flow  rates  and 
physical  dimensions.  Coupled  with  a  software  that  optimizes  cost,  the  designer  can  specify  the 
most  cost  and  performance  efficient  plant,  or  best  value  for  the  customer. 
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N/N0  (LOG  Inactivation  Achieved) 


Temperature  D 


Conclusion: 


There  is  a  need  for  a  method  of  predicting  and  optimizing  inactivation  of  pathogens  in 
drinking  water  treatment  plants.  In  the  interest  of  public  health,  new,  more  stringent  standards 
must  be  implemented,  achieved,  and  enforced.  The  design  and  operation  of  a  bubble-diffuser 
ozone  contactor  should  not  be  strictly  based  on  comfort,  cost  or  tradition.  As  more  is  discovered 
about  Cryptosporidium  spp.,  disinfection  and  filtration  should  jointly  and  intelligently  be  utilized 
to  achieve  the  99.9999%  inactivation  level  desired.  The  model  described  in  this  paper  is  based 
on  mass  balance  principles  and  has  been  validated  by  experimental  results.  When  incorporated 
in  a  spreadsheet  type  software,  the  equations  take  the  form  of  a  simple,  user-friendly  application 
with  graphically  displayed  results.  It  provides  predictions  of  dissolved  ozone  concentrations  and 
the  corresponding  pathogen  inactivation  which  could  be  reasonably  anticipated  for  a  given  set  of 
physical  characteristics  and  operating  parameters  with  a  very  high  level  of  accuracy.  The 
application  to  design  and  operation  of  bubble-diffuser  ozone  contactors  is  unlimited  and  at  a 
minimum  provides  excellent  guidance  and  recommendations  in  order  to  achieve  optimum 


inactivation. 
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Appendix 


Appendix  I.  Notation  and  General  Calculations 


Symbol  Definition 

A  Chamber  Area  (cm2) 

a  Effective  Area  (cm1) 

Cl  Concentration  of  ozone  in 

liquid  phase  (mg/L) 

Cl*  Interfacial  ozone  concentration  (mg/L) 

Cg  Concentration  of  ozone  in 

the  gas  phase  (mg/L) 

D  Chamber  Diameter  (cm) 


Formula 


rix 


D 


6xUr. 


Wb  -  Ul) 


d 

Dispersion  Number 

dB 

Average  Gas  Bubble 
Diameter  (cm) 

dBO 

Gas  Bubble  Diameter 
at  Diffuser  Outlet  (cm) 

Eg 

Turbulent  Gas 
Diffusivity  (cm2/sec) 

El 

Turbulent  Liquid 
Diffusivity  (cm2/sec) 

El 

Ul-L 

d bo  0-2 1 UG 


6E  -8x- 


Gi„  Percent  of  Ozone  in  Applied 

Gas  by  Weight 


Kd  Ozone  Decomposition  Rate  (sec'1)  1 .808£14  x  Exp\ 


Kl 


Kol 


Kn 


L 

m 

M 

N 


Nd 


Nl 


Nn 


Sh*PL 

Liquid  Mass  Transfer  d 

Coefficient  (cm/sec) 

Overall  Mass  Transfer 
Coefficient  (cm/sec) 

Inactivation  Constant  for 
Pathogen  (mg/L-min) 

Chamber  Height  (m) 

Henry’s  Constant  Log  (m)  =  3.25 

Mass  (mg) 


Quantity  of  Oocysts 

Dimensionless  Ozone 
Decomposition  Rate 

Dimensionless  Mass  Transfer 
Coefficient 


kd-l 

UL 


KL  x ax L 

UL 


Dimensionless  Inactivation  Constant 
for  Pathogen 


Qg  Gas  Flow  Rate  (L/min) 

Ql  Liquid  Flow  Rate  (L/min) 


-9746 > 
286.15, 


840 

T 
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rc  Mass  reaction  rate  (mg/L) 


fL 


Re 

Sc 


Sh 


s 

t 

T 

T 


UG 


UL 


M-l 


Vb 


V 

X 


Mass  transfer  rate  (mg/L) 
Reyolds  Number 

Schmidt  Number 

Shaffer  Number 

Stripping  Factor 
Time  (sec) 

Temperature  (°K) 


Vg  X  <jjg 
VL 


vL 


Dl 


2 +  .0187 


Rem 


x,S; 


^.072  “ 

D667 

\  UL  ) 

m-Qg 

Ql 


t  + 273.15 


Temperature  (°F) 

Gas  Flow  Velocity  (cm/sec) 

Liquid  Flow  Velocity  (cm/sec) 
Dynamic  Viscosity  (centipoise) 

Bubble  Rise  Velocity  (cm/sec) 


Qg_ 

A 

Ql 

A 

_ 100 _ 

t  -  8.4385+  ^80784  +  (/  -  8.43  5  J  \ - 1 20 


Kinematic  Viscosity  (cm2/sec) 
Length  (cm) 


Ml. 

Pl 


44 


Z  Dimensionless  Length  Coefficient 

A  Denotes  “change  in”  the 

associated  variable 

d  Denotes  “infinitesimal  change  in” 

the  associated  variable 

s  Liquid  Volume  Fraction 

0  Mass  Flux 

0 


Dimensionless  Time  Coefficient 


Appendix  II.  Mass  Balance  for  Counter-Current  Ozone 
Chamber 


No 


Liquid  Side: 

Mass  Accumulation  =  Mass  In  -  Mass  Out  +  Mass  Reacted  +  Mass  Transferred 
A M  -  O  -  (d>  +  AO)+  V  -rL+V -rc 


AM  =  Cl -A-s-AX 


9.-Bl.A.s.%  +  Ql.Cl 


rL=KOL-a-e-(cl-CL) 

rc  =  —^D  ■  £ '  CL 
So  the  equation  becomes: 


A(AA/)  =  -AO- At  +  V  -rL  -A t  +  V ■ rc  -At 

(  A 2C  ^  /  \ 

ACLAs-AX  =  ElAs - l--QlACl  At  +  As-AX-KOLa(C*-CL)At-As-AX-KL 

l  AX  , 


Dividing  through  by  A-s-AX -At  results  in: 


CLAt 
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Gas  Side: 


Mass  Accumulation  =  Mass  In  -  Mass  Out  +  Mass  Reacted  +  Mass  Transferred 
AM  =  Cg -A-(l-e)-AX 

®  =  -Ea-A(l-e)~£--Qc-Cc 

rL  =  -a-e-(cL  —  CL) 

rc=-KD-fy-s)CL 
So  the  equation  becomes: 

A(AM)  =  -AO  •  At  +  V  ■  rL  •  At  +  V  ■  rc  •  At 


ACgA{\-s)-AX=\  EG4s)^+aACG  \st+A(\-£)-&X-KLc(CL  -C,  )&-A{l-£)- AX- KDCGAt 
\  A^r 


Dividing  through  by  A  ■  (l  -  s )  •  AX  •  At  results  in: 


acg  a2cg  t  Qg  A Cg _ s_ 

At  G‘  AX2  A-(l-s)  AX  (l-f) 


KL-a-(ci-CL)-KD-C0 


Taking  the  limit  as  the  delta  becomes  infinitely  small,  A  ->•  0  results  in: 


dCr  „  d2Cr  Ur.  3C, 


'G  _  j?  _  _ _ 

a  G  ’  ax2  "i-s ’  ax  (l-s) 


■KL-a{cl-CL)-KD-Cc 


To  non-dimensionalize  the  equation,  multiply  each  term  by  —  : 

U  j 


L  dCr.  En  ■  L  dlCn  Ug-L  dCG  K,-a-L-e 


_^jjL±±.(f.cyKjLL 

/I  TT  V  L  L /  JJ  0 


uL  a  uL  ax2  (1  -e)-uL  ax  fy-e)-uL 

The  non-dimensionalized  equation  then  becomes: 


dCG  _  Eg  d  Cq 


Ur. 


dCG  Nr  s 


0©  Ul-L  dZ 2  UL-(l-s)  dZ  (l-e) 


-ttA-{c‘l-cl)-ndcc 
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3C 

If  the  contactor  is  allowed  to  reach  steady-state,  then  — —  =  0 . 

50 

Additionally,  considering  that  diffusivity  in  the  gas  bubbles  is  extremely  small  and  ozone 
is  non-reactive  in  air,  we  can  eliminate  two  additional  terms. 

The  mass  balance  for  the  gas  side  becomes: 


UG  dCG  Nr  •£ 


UL-(l-s)  dZ  (1-s) 

C  Wl  *  U 

Since  C*L  =  — ,  and  S  = - — ,  the  equation  becomes, 


U, 


fC^ 


V  m  j 


(l-ff)  5Z  (l-s) 


Nr 


£o_c 

'-L 

m 


=  0 


Now  we  have  two  equations  representing  our  system.  Assuming  e  =  1  for  ozone,  the  two 
equations  reduce  to  the  final  form: 

WCl-o 


dZ  dZ 


m 


II. 


fC  ^ 


V  m  J 


-MSl-c. 


=  0 


dZ  S  y  m  j 
Using  the  boundary  conditions: 
Liquid: 

CI(Z=0)  =  CUJ,+rf 


dZ 


lz=o 


(Inlet) 


5Q 

az 


^  =  0  (Outlet) 


Gas: 


C0(Z-1)=C0, 


IN 
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If  in  Equation  II,  we  let  Y  =  — ,  we  get: 

m 


dZ  S'  L’ 

We  can  rearrange  and  solve  for  Cl  in  terms  of  Y  as  follows: 

CL  =Y — — 

L  Nl  dZ 

With  an  expression  for  Cl  in  terms  of  Y,  we  can  develop  expressions  for  the  first  and 
second  derivatives  of  Cl  as  follows: 


ecL  _  dY  s  a2r 

dZ  ~  dZ  Nl  '  dZ2 


and 


d2CL  _d2Y  _  S  tfY_ 
dZ2  ~  dZ2  Nl  '  dZ3 

Substitution  into  Equation  I  results  in: 


f a2Y  s  a3r^i 


8Z2  Nl  dZ 
or 


(dY  s  d2r^ 
dZ  Nl  '  dZ2 


■ 

+  *i 

Y- 

) 

V 

S  dY 
NL'dZ 


-Nr 


ry  S  dY^ 
N  L  '  dZ 


=  0 


d2Y 


.  S  d3Y  dY  S  d2Y  _  dY  ..  S  dY_n 

dz2  N,  dZ3  dZ  Nl  dZ2  dZ  Nl  dZ 

Combining  like  terms  results  in: 


d-S 

d3Y 
*  -  + 

f,  s ) 

d2Y 

NL 

dZ3 

l  Nl) 

dZ2 

l  Nl  ) 

8Z  ° 


We  now  have  a  linear  differential  equation  with  constant  coefficients. 
Using  a  particular  solution  of  Y  =  EXP(AZ),  we  get: 


d-S 


A3  • EXP(AZ)+ 


d  +  - 


N 


A2  • EXP{AZ)+ 


L  ) 


N°  ■ S  - 1  -  s  }a  •  exp{az )  -  nd  •  exp(az) 


V 


NL 
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The  general  solution  to  this  equation  is: 

Y  =  a,  •  EXP(A]Z)+cc2  •  EXP(A2Z)+a3  •  EXP(A3Z) 


dY  d2Y  d3Y 

With  this  solution,  we  can  develop  expressions  for  — , — jand  — -  which  will  be 

oZ  dZ  oZ 

useful  later  in  the  incorporation  of  boundary  conditions. 

—  =  a,  •  JLi  •  EXP(AlZ)+ a2 -A2  ■  EXP(A2Z)+  a3  •  A3  ■  EXP(A3Z ) 
dZ 


£Y 

3Z2 


=  arA 2  ■  EXP(A1Z)+  a2  •  A2  •  EXP(A2Z)+  a3  •  A]  ■  EXP(A3Z) 


d_Y_ 
dZ 3 


=  arA ( •  EXP{AxZ)+  a2  ■  A\  ■  EXP{A2Z)+  cc3-A\  •  EXP(A3Z ) 


XT  ,  J  ,  .  .  _  dCL  d2CL 

Now  we  can  also  develop  expressions  tor  CL, - and - — : 

dZ  dZ 


NLdZ 


Q  ” 


EXF{AlZ)-jj-AlEXl{AxZ) 


a,  + 


EXF(A2Z)-±-A2EXI{A2Z) 


kXi  + 


exf{a3z)-^-a3exp{a3z) 


ax 


5Q 

dZ 


AEXFiAZ)  -  %EXF{\Z)  [or,  + 


A2EXF{A2Z)  - — A^EXfiAZ) 

Nl 


Aexf{az)  -  jj-  Aexf{az) 


d2CL 
dZ 2 


Aexp(az)-—%exp(az) 


N, 


x\exp{a2z)~ — a\exp(a2z) 


N, 


«,  + 


Aexp(az)-—a]exp(az) 


N, 


a , 


We  can  now  apply  the  first  boundary  condition: 

0^=0)=^  L 
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f  s  '  ' 


1-— 4 

V  J 


Fl  + 


1-— A, 

V  J 


( 


F2  + 


V  *LJ 


F»  = 


M  dS  2 
“A  “TT^i 
V  ^  / 


K  + 


c^-- — 4  k  + 


L  J 


<^3  ,,  A  °3 

■/'£  / 


Rearranging  results  in: 
Equation  1 : 


1-— X[-dXi+  — 


CTj  + 


^5 


V  *L 


Nl 


dS  ' 


1  -Z-^-d^+^-%  a2  l-ZL.A3-dA3+^ 


W, 


a3=0 


The  second  boundary  condition  is: 


dC, 

dz  12 

Equation  2 


=  =  0  resulting  in: 


lEXfay-Z-tEXfa) 


N, 


«.+ 


A'Enfa)-±%EXfa) 


N, 


02  + 


XJEXfa)-±%EXfy ,) 


^3=0 


For  the  third  boundary  condition,  we  need  an  expression  for  Cg  in  terms  of  CL.  From 
Equation  I,  we  can  solve  for  Cg  as  follows: 


Q=- 


^  ^  fND-m  ' 


d-m  d  CL  m  dCL 


Nl  dZ2  Nl  dZ 


+ 


nl 


■  +  m 


CL 


Now  we  can  apply  the  boundary  condition  :  CG  (Z  =  l)  =  Q?,  IN 


m 

~N~, 


| 'x]EXP  (/l,)-  -jj-A\EXP  (A,)ja,  +  j  X\EXP  (A2)-  -jL-X \EXP  (A2)ja2  +  ^EXP  (A,)- -^-A’£XP  (A3)ja3 

\  /  o  \  f  C 

XlEXP(Al)-—A]EXP(Xi)  a,  +  A1EXP(A2)-—A21EXP(A1)\a1  +  A,EXP  (A3)~  —  A\EXP  (As) 

"L  J  V  Nl  )  \  Nl 


Npm 

Nr 


+  m 


EXP  (A, )  -  A.  EXP  (A,  )la,  +  {eXP{A1)-^-A1EXP{A1)^1  +  f  EXP(A3)--^-A3EXP(A3) 

nl  J  f  nl  j  y  nl 


=  Cr. 
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Combining  like  terms  results  in: 
Equation  3: 


d’tn  (  ,2  rvD^  \  S  -3  T 7Vr>(«  \ 1  .  m  (  1  Z7  VD  ( 2  \  &  rw{  1  \  I  -L  (  N  pm_  ¥  t?YL>(  3  \ _ ^ 


X\EXP(Xx)-^-X\EXP{Xx)U^-  XxEXP(Xx)-^-X]EXP{Xx) 


+  — +  m 


EXP  (Xx )  -  X,  EXP  (Xx )  a,  + 


X]EXP(X2)-4-X\EXP(X2)  +  -£-  X2EXP(X2)-4-X22EXP(X2)\  +  \  ^-+m  EXP(X2)--±-X2EXP(X2)  a 


-^L\x]EXP(X2)--^-XlEXP(X2)  +-ZL  X2EXP(X3)--^-X\EXP(X3)\  +  \  —2—+m  EXP(X3)--£-X}EXP(X2)  a2=C, 


We  now  have  three  equations  with  three  unknowns:  ax  ,a2  and  a3 .  (Once  we  solve  for 
the  roots  A,  ,  A2  and  A3 .) 

Recalling  that  the  roots  can  be  determined  from  the  linear  differential  equation: 


d  S ^  A3+U  +  — M24^^-l-S|-A-jVn  =0 


With  the  roots  A, ,  A2  and  A3 ,  our  three  equations  can  be  solved  simultaneously  for  ax ,  a2 


and  a3 . 


S  3Y 

Recalling  that  C,  =Y - ,  we  can  solve  for  the  concentration  of  dissolved  ozone 

L  Nl  8Z 

at  any  point  in  the  contactor  with  the  following  equation: 


CL(z)=  EXP  (Xxz)-  —XxEXP  (XxZ)  la,  +  EXP  (X2Z)~  -f-X2EXP  (X2Z)  \a2  +  EXP  (X3Z)-  —X3EXP  (X3Z)\a 
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Appendix  III.  Mass  Balance  for  Co-Current  Ozone  Chamber 


Liquid  Side: 

Mass  Accumulation  =  Mass  In  -  Mass  Out  +  Mass  Reacted  +  Mass  Transferred 
A M  =  ®  -  (®  +  A®)+  V  ■  rL  +  V  ■  rc 


A M  =  Cl -A-s-AX 

rL  ~  K0L  '  a  '  S  '  {pL  “  CL  ) 

f'c  =  ^ D  '  ^  ^ L 

So  the  equation  becomes: 


A(A M)  =  -  A®  •  At  +  V  -  rL  •  At  +  V  ■  rc  -  At 


ACLAs-AX 


ElAs 


A 2Q 

AX 


-Ql^l 


At  +  As  ■  AX  ■  KOLa{c*L  -  CL  )At  -  As  -  AX  ■  KDCLAt 

J 


Dividing  through  by  A-s-AX -At  results  in: 
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A Q 
At 


=  EL-^j--^-^  +  KOL-a-(CL-CL)-KD-CL 
AX2  A-s  AX  v  > 


Taking  the  limit  as  the  delta  becomes  infinitely  small,  A  -»  0,  and  additionally 
substituting  Kl  for  Kol  since  the  mass  transfer  resistance  on  the  gas  side  is  negligible 
results  in: 


X'-=E,-^-^-^  +  KL-a-(cl-CL)-KD-CL 


dt 


dX1  s  dX 


To  non-dimensionalize  the  equation,  multiply  each  term  by 


UL 


±_  SCL  El  •  L  cEC_l  UrL  dCL  ,  KL-a-L  /  .  x  KD-L 

UL'  a  UL  '  dX2  e-UL  dX  u,  \  L  L)  T1  L 


u, 


Non-dimensionalized  coefficients  will  be  defied  as  follows: 


0  = 


t-uL 


Kl  • a-L 
UL 


Z  =  * 
L 


N 

D  U, 


Ul-L 


r 

r*  - 
W  ~ 

m 


The  non-dimensionalized  equation  then  becomes: 


dCL  =d.^_l.^x+NL\£o_.cL  I -nd-cl 


50 


dZ2  s  dZ 


m 


dC 

If  the  contactor  is  allowed  to  reach  steady-state,  then  — —  =  0  and: 

50 


The  mass  balance  for  the  liquid  side  becomes 

fCt 


,|  ^2._  cL 


dZ 2  s  dZ 


m 


-Nd-Cl=  0 
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Gas  Side: 


Mass  Accumulation  =  Mass  In  -  Mass  Out  +  Mass  Reacted  +  Mass  Transferred 
AM  =  CG -A-(l-s)AX 

«  =  -EGA4-e).^-  +  QG-Ca 

rL=KL-a-e-(c*L-CL) 

rc  ~  •  (l  —  £■)  •  Cl 

So  the  equation  becomes: 

A(A M )  =  -AO  •  At  -  V  •  rL  ■  At  +  V  •  rc  •  At 


ac04i-£)-ay=  EG4 W-4i-£)-&r -kl<(cl  -c>-4i-s> at  KDcG6t 

V  A X  J 


Dividing  through  by  A  AX  -At  results  in: 


ACg  =E  A2Cg  Qg  ACq _ £ 


-KL-a-{c:-CL)-KD-Cc 


At  G  AX 2  A-{l-e)  AX  (l-s) 

Taking  the  limit  as  the  delta  becomes  infinitely  small,  A  0  results  in: 


£k  =  jE  UG  dCG _ e_ 

a  G  ax2  1  -s  ax  (i-f) 


■KL-a(c[-CL)-KD-CG 


To  non-dimensionalize  the  equation,  multiply  each  term  by 


UT 


L  dCr  Er-L  d1Cr  U r  -L  SCG  K,  -a-L-s 


_ _ _ _ —Xc*  -C  )-  —  —  •  c 

uL  a  uL  ax1  (i -s)-uL  ax  (1  -s)-uL  ^  L  1  uL  G 

The  non-dimensionalized  equation  then  becomes: 


5CV.  Er  d2Cr, 


Ur. 


dCr  N,  s 


a©  ul-l  az2  uL(i-e)  az  (i-^) 


(c;-q)-jvQ 


56 


If  the  contactor  is  allowed  to  reach  steady-state,  then 


50 


=  0. 


Additionally,  considering  that  diffusivity  in  the  gas  bubbles  is  extremely  small  and  ozone 
is  non-reactive  in  air,  we  can  eliminate  two  additional  terms. 

The  mass  balance  for  the  gas  side  becomes: 


Ur  5Cr.  N,  s 


UL  •(!-*)  dZ  (l-a) 


^•(c;-ct)=o 


Since  C[  =  — ,  and  S  =  —  — 
m  UL 


,  the  equation  becomes, 


r  n  \ 


V  m  J 


7 - x - +  T^—\‘Nl  ’[—-Cl 

(l -s)  dZ  (l-e)  V  ”1  J 


=  0 


Now  we  have  two  equations  representing  our  system.  Assuming  e  =  1  for  ozone,  the  two 
equations  reduce  to  the  final  form: 


I. 


dCL 
dZ 2  dZ 


+NL 


'—-cl)-nd-cl=0 


m 


II. 


(r  \ 


)  xnl 

dZ  S 


~Cr 


=  0 


J 

Using  the  boundary  conditions: 
Liquid: 


C i  (Z  —  0)  -  CLIN  +  d  |z=0  (Inlet) 


5Q  | 

dZ  lz 


=)  =  0  (Outlet) 


Gas: 


Co(Z  =  0)=Co,„ 
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If  in  Equation  II,  we  let  Y  =  — — ,  we  get: 

m 


— +— (f -q)=o 

ez  s  ' 

We  can  rearrange  and  solve  for  Cl  in  terms  of  Y  as  follows: 

^  v  S  8Y 

Ci  —  Y-\ - 

1  Nl  8Z 

With  an  expression  for  Cl  in  terms  of  Y,  we  can  develop  expressions  for  the  first  and 
second  derivatives  of  Cl  as  follows: 


8C,  dY  S  d2Y 

— ±  =  —  + - -  and 

dZ  8Z  N,  8Z2 


82CL  =  d2Y  S  &Y 
8Z2  ~  8Z2  +  Nl  '  8Z 3 


Substitution  into  Equation  I  results  in: 


d 


f  82Y  S 

a3r^| 

f  8Y 

S 

J _ 

82Y 

j 

laz2  nl 

az3J 

[dZ 

Nl 

'  8Z2 

or 

d2Y+d.  S 

82Y 

8Y 

S 

82Y 

8Z2  N, 

8Z 3 

8Z 

NL 

8Z2 

\ 

f  5 

8YS 

r 

+  nl 
J 

Y- 

Y  + - 

l  nl 

8Z  J 

~  N  d 

\ 

Y  + 


S_  8Y 
NL'dZ 


=  0 


ar 

8Z 


S  8Y_ 
NL‘  8Z 


Combining  like  terms  results  in: 


d-S 

83Y 
- r  + 

f,-Al 

82Y  r 

Nl 

az3 

«l) 

'  az2  [ 

Nn-S 


1  +  S  +  — 


N 


L  J 


8Y 

8Z 


-Nd-Y  =  0 


We  now  have  a  linear  differential  equation  with  constant  coefficients. 
Using  a  particular  solution  of  Y  =  EXP(XZ ),  we  get: 

S 


'd-S' 


A3  •  EXP(XZ)+ 


d-- 


N, 


U 2  -EXP{XZ)~ 


f  N  -S^ 
1 +  S  +  ^ — 
NL 


a  ■  exp(xz)-nd  ■  exp(xz) 
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The  general  solution  to  this  equation  is: 

Y  =  ar  EXP(AlZ)+  a2  ■  EXP{A2Z)+a3  ■  EXP{A3Z) 


BY  B2Y  83Y 

With  this  solution,  we  can  develop  expressions  for  — , — -and  — -  which  will  be 

BZ  BZ  BZ 

useful  later  in  the  incorporation  of  boundary  conditions. 

—  =  a,  •  V  EXP(AlZ)+  cc2-A2-  EXP(A2Z)+a3  ■  A3  ■  EXP(A3Z) 
dZ 


dx 

ez2 


=  a,  -A2  ■  EXP(A1Z)+  a2  ■  A\  •  EXP(A2Z)+  a3-A2- EXP(A3Z ) 


BX 

BZ3 


=  ax-A\-  EXP{A]Z)+  a2  ■  A\  ■  EXP{A2Z)+  cc3-A\ •  EXP(A3Z) 


.  r  „  BC,  ,  82C, 

Now  we  can  also  develop  expressions  for  CL ,  — —  and  : 

BZ  BZ 


Cl=Y  + 


S  BY 

nl  bz 


Q 


EXl{AlZ)+—A[EXl{AlZ)  |aj  + 


EXl{A2Z) + -^-A1EXl{A2Z) 


a,  + 


EXl{A3Z) + j^A3EXl{AiZ) 


Vx3 


bcl= 

BZ 


AiEXfyxZ)+jjr  %EXfaz) 


Fi+\ 


a2exi(a2z)+— fiExfaz) 


N, 


M 


d2Ch 

BZ2 


%EX$AlZ)+—  %EXfaz) 


Nr 


k  + 


X1ex^a2z)+—^ex^a2z) 

N i . 


A \EXfatil+±-%EXfa$ 


F3 


We  can  now  apply  the  first  boundary  condition: 


cI(z  =  o)»</^- 


Z= 0 
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(  s  '  t 


1+— 4 

.  *LJ 


Fi  + 


V  ^ 


6^2  + 


i+^ 


V 


.,  ,  dS  ,  f  ,  dS  , 
<A+— ^1  «,+  ^2  +TT^2 
Nl  )  \  Nl 


«2  + 


J^3+— ^ 


a, 


L  J 


Rearranging  results  in: 
Equation  1 : 


f  c  ,/c 

lH - Z,  —d/ lj - ^  la,  + 

V  Nl  ‘J 


A 


.  5  .  dS  2 

1  H - >^2  ~  UA^ - ^2 

V  J 


^  C  .70  ^ 

1  +— A3-dZ3-  — Z]  a3=0 
v  ^  J 


The  second  boundary  condition  is: 


5Q 

az 


=  0  resulting  in: 


Equation  2: 


la,  + 


Z2EXl{Z2)+  — ■%EXfa)a2  + 


N, 


Z3EXl{Z3)+-r^EXl{Z3) 


Nr 


\a3  =0 


For  the  third  boundary  condition,  we  need  an  expression  for  Cg  in  terms  of  Cl-  From 
Equation  I,  we  can  solve  for  Cg  as  follows: 


d-m  d2CL  m  dCL  f ND-m  ^ 

G  nl  dz2  nl  dz  [  nl  j 


■Q 


Now  we  can  apply  the  boundary  condition  :  CG  (Z  =  0)  =  CGIN 


d  *  m 


NL 


'a* EXP  (A, )+  —  A] EXP  (A, )1«,  +  f A\EXP  (A2 )+  -f-  A\EXP  (A2 )la2  +  [  A*  EXP  (A, )+  A\EXP  (A, ) 
NL  J  \  N  l  )  \  nl 


m 

A ~l 


AtEXP  (A,  )+—AfEXP(Al))a,+(Ai  EXP  (A2)+^-A22  EXP  (A2)k+f  A3  EXP  (l, )+  A*  EXP  (l3 ) 
Aft  J  ^  J  ^  > 


V  ^ 


+  m 


\  / 

a,  + 

y  v 


£Ay(21)+-^A1£XP(ll)|a,+|  £XP(A2)+-|-A2£ZP(^2))«2+|  EXP{A2)+±-A2EXP(Ai) 


N, 


=  Cn 
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Combining  like  terms  results  in: 
Equation  3: 


N, 


d -m 

'~n7 


X^EXP(Xl)+  4-^EXPm)  +  M  X, EXP&) +  -§-%EXPM )  +  fe  +  mY  £H>(A,)+  -g-A,EXPM 

nl  J  nl\  nl  J  \  ni  A  ^ 

(  %EXP(Z2)+  ^-X\EXP(X2) ]  +  ■%-(  X2EXP(X2)+  JL. x\exp{x2 ) 

V  Ni  )  NiV  nl 


s^+m  \  EXP(X2)+—X2EXP(X2) 
Ml  j 


a,  + 


a,  + 


dm 

"aT 


A/*  J  ArL^  A/£  J  v  A  ^  /_ 


We  now  have  three  equations  with  three  unknowns:  ax,a2  and  a3 .  (Once  we  solve  for 
the  roots  /l,  ,/l2  and  A3 .) 

Recalling  that  the  roots  can  be  determined  from  the  linear  differential  equation: 


f d-S A 


Nl 


■A3  + 


d  — 


N 


f 


■A2  - 


LJ 


1  +  5  + 


V 


No± 

Nr 


■a-nd=  0 


With  the  roots  AX,A2  and  A3 ,  our  three  equations  can  be  solved  simultaneously  for  ax,a2 
and  a3 . 


S  dY 


Recalling  that  Cr  =  7  + - ,  we  can  solve  for  the  concentration  of  dissolved  ozone 

5  L  Nl  dZ 

at  any  point  in  the  contactor  with  the  following  equation: 


(z)=  [  EXP  (A,z)+  —  X,EXP  (t.Z)]**,  +  { EXP  (X2z)+ X2EXP  (t2Z)L  +  [ EXP  (X3z)+  X,EXP  (l, 

V  M L  )  \  Nl  J  y  Jy  L 
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Appendix  IV.  Mass  Balance  for  Reactive  Ozone  Chamber 


No 


Mass  Accumulation  =  Mass  In  -  Mass  Out  +  Mass  Reacted  +  Mass  Transferred 
A M  =  O  -  (d>  +  A®)+  V-rL+V  -rc 


AM  -CL -A-s-AX 

^ :  ~El  '  A  -  s  -  +  QL  -  CL 

rL  =  0  (No  ozone  vapor  bubbles) 
rc  =  —Kd  •  s-Cl 
So  the  equation  becomes: 

A(AM )  =  -A<I>  -At  +  V-rc -At 


ACLAs-AX  = 


AC 

ELAs-^--QLACL 


to- As -  AX  -  KDCLAt 


Dividing  through  by  A-s-AX -At  results  in: 
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A Q 

At 


A 2C 


—  p - L  Z-'L - L 

L  A  \r2  A  ~  A  V 


AZ 


AC, 


AZ 


-Z„-C, 


Taking  the  limit  as  the  delta  becomes  infinitely  small,  A  -»  0 ,  and  additionally 
substituting  Kl  for  Kol  since  the  mass  transfer  resistance  on  the  gas  side  is  negligible 
results  in: 


SC±  =  E  <?2Q 

a  L  ax1 


kd.q 
e  ax  D  1 


To  non-dimensionalize  the  equation,  multiply  each  term  by  — —  : 


_L_  SC \  <?2Q  U^L  dCL  KpJL 

JFL  ’  a  ~  uL  '  ax2  s-uL  ax  uL  L 

Non-dimensionalized  coefficients  will  be  defied  as  follows: 


©  = 


t-uL 


L 


°  u, 


d. 


Ul-L 


The  non-dimensionalized  equation  then  becomes: 


.>£l.Nb.Ci 

5© 


dZ 2  e  dZ 


8C 

If  the  contactor  is  allowed  to  reach  steady-state,  then  — -f-  =  0  and: 

50 

The  mass  balance  for  the  liquid  side  becomes: 


d 


d2C,  1  dC, 


-Nd-Cl=  0 


dZ2  s  dZ 

Again,  assuming  £  =  1  for  ozone,  the  equation  reduces  to  the  final  form: 
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l-nd-cl=  0 


J  d2C,  dC 
dZ 2  dZ 

If  we  let  Y=Cl,  we  get: 


d 


a2r  8Y 

dZ2  dZ 


-ndy  =  o 


We  again  have  a  linear  differential  equation  with  constant  coefficients.  Using  a  particular 
solution  of  Y  =  EXP(A  •  Z),  we  get: 


d  ■  A2EXP(A  ■  Z)  -  AEXP(A  •  Z)  -  ndexp(a  •  Z)  =  0 
After  dividing  through  by  EXP(A  ■  Z),  we  get: 


d-A2 -A-ND  =0 

The  general  solution  to  this  equation  is: 


CL=Y  =  f3r EXP(A4  • Z)+p2  •  EXP(A5  •  Z) 


We  can  now  develop  an  expression  for  —  which  will  be  useful  later  in  the  evaluation  of 

dZ 

boundary  conditions: 

~ ■=  fi,  -K  ■  exp(K •  z)+ A-V •  z) 

oZ  dZ 

Our  first  boundary  condition  is  that  the  input  to  the  reactive  chamber  is  the  output  of  the 
previous  chamber: 

CL  (Z  =  0)  =  (Inlet)  Since  CL=Y,  we  get: 


Equation  1 :  Px+P2=  Cv~chamber 

The  second  boundary  condition  is: 


dC, 

dZ 


Zml  =0,  (Outlet) 


resulting  in: 
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Equation  2:  \A4  •  EXP(Aa  )]/?,  +  [As  ■  EXP(A5  )]/?2  =  0 


We  now  have  two  equations  with  two  unknowns:  /?,  and  (52 .  (Once  we  solve  for  the 
roots  A4  and  As .) 

Recalling  that  the  roots  can  be  determined  from  the  linear  differential  equation: 
d-A2  -A-Nd=0 

With  the  roots  A4  and  As ,  our  two  equations  can  be  solved  simultaneously  for  /?,  and 

Pi- 

We  can  now  solve  for  the  concentration  of  dissolved  ozone  at  any  point  in  the  reactive 
chamber: 


CL  (Z)=  px  ■  EXP  (A4  ■  Z)+  J32  ■  EXP  (As  •  Z) 


Since  no  gas  bubbles  are  applied  in  a  reactive  chamber,  it  may  be  assumed  that  turbulent 
dispersivity  is  minimal  and  approaches  “plug  flow”  conditions.  Thus,  the  dispersion  term 
may  be  ommited  and  our  equations  simplified. 

The  mass  balance  for  the  liquid  side  becomes: 


dC, 

dZ 


+  Nd 


•Cl 


=  0 


Using  a  particular  solution  of: 

CL  =y  =  exp(a-z) 

We  get  the  following  equation: 

2  •  EXP(A  -Z)+Nd-  EXP{A  ■  Z)  =  0 


Solving,  we  find  that  the  root  ( A )  =  -No- 


The  general  solution  to  this  equation  is  then: 
CL  =Y  =  cc-EXP(-Nd-Z ) 
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Only  one  boundary  condition  is  required.  The  input  to  the  reactive  chamber  is  the  output 
of  the  previous  chamber: 


CL  (z  =  o)  =  c““  (Inlet) 


This  allows  for  the  development  of  an  equation  representing  the  concentration  of 
dissolved  ozone  at  any  point  in  the  reactive  chamber  when  dispersion  is  assumed  to  be 
negligible: 

CL  (Z)  =  •  EXP(~  nd  •  z) 
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Appendix  V.  Mass  Balance  for  Cryptosporidium 
Inactivation 


No 


Ax 


0 

0+ A© 


Qs 


N 


Mass  Accumulation  =  Mass  In  -  Mass  Out  +  Mass  Reacted 
AM  =  O  -  (O  +  AO)+  V  •  rN 

AM  =  N  ■  A-  AX  (Where  N  =  Quantity  of  active  oocysts) 

<f--EL-A~  +  QL-N 

rc=-K„-N-CL 
So  the  equation  becomes: 


A(A M )  =  -AO  •  At  +  V  •  rN  ■  At 


AN -A- AX  = 


\ 


ElA--QlAN 
1  NX 


At-A-AX-K„-N-CrAt 


Dividing  through  by  A  •  AX  ■  At  results  in: 


AN  „  A2N  Ql  AN 


At 


=  E,  ■  =L^r  -  —  •  —  -Kn-N  -CL 
L  AX2  A  AX  N  1 


Taking  the  limit  as  the  delta  becomes  infinitely  small,  A  -»  0 ,  results  in: 


—  =  El -^-Ul  ~-Kn  ■ n-cl 

a  ax2  ax 


To  non-dimensionalize  the  equation,  multiply  each  term  by 


Ur 


_L_  aN  _El-L  a2N  Uj^L  cN_  Q 

uL'  a  ~  uL  '  ax2  uL  '  ax  uL  L 

Non-dimensionalized  coefficients  will  be  defied  as  follows: 


^  t-u, 

©  = - — 


Ul-L 


z  =  * 

L 


N  =K"'L 

JV  N 


Ur 


The  non-dimensionalized  equation  then  becomes: 


80  dZ2  8Z 


dN 

If  the  contactor  is  allowed  to  reach  steady-state,  then  —  =  0  and  in  final  form,  we  have: 

5© 


d- 


82N  8N 


- Nn  ■  N -CL=0 


8Z2  8Z 
Our  boundary  conditions  are: 


N(Z=0)  =  N0  +  d 


a_  N 
a  z 


I  2=0 


(Inlet)  and 


8N_ 

8Z 


lz=l 


=  0 


(Outlet) 
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Using  a  finite  difference  approach,  we  can  develop  a  profile  for  Cryptosporidium 
inactivation  through  the  contactor. 


Z=  1 


-m- 


iz=  1 
iz  =  2 
iz  =  3 
iz  =  4 
iz  =  5 


iz  =  w-3 
iz  =  w-2 
iz  =  w-1 
iz  =  w 


iz  -  1 


W 


I 

Space  increment  A  z  =  — j- 


If  we  let  w  =  401,  A  z  =  0.0025 


^  JV(iz  +  l)-W(/z-l) 

The  expression  for  — •[,  = - — - 

At  the  top  of  the  contactor  (where  iz  =  1) ,  we  have: 

g  N .  _  -iV(/z  +  2)  +  4 N{iz  + 1)  -  3N(iz) 

~g~Z'iz=l  2AZ 

At  the  bottom  of  the  contactor  (where  iz  =  w ) ,  we  have: 


g  N  |  _  W(fe-2)-4AT(fe-l)  +  3tf(fc) 

2AZ 

Ou  expression  for  every  increment  in  between  is: 

g2N  _  + 1)  -  2N(iz)  +  N(iz  - 1) 

~d  Z2 (AZ)2 

At  iz  =  1  and  iz  =  w ,  we  can  apply  our  boundary  conditions: 
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At  iz  =  1 : 


N(l)-d 


-N(3)+4N(2)-3N(i)' 

2AZ 


=  N„ 


At  iz  =  w: 

Njw  -  2)  -  4N(w  - 1) + 3A(w)  _ 
2AZ 


=  0 


By  rearranging  ,  the  first  equation  becomes: 

After  grouping  like  terms,  the  first  equation  becomes: 


1) 


1  + 


3d 

2AZ 


N(l)-^N(2)+^=N(3)^N, 
AZ  2AZ 


B(l)  D(l)  X(l) 

Similarly,  the  second  equation  becomes: 


G(l) 


2)  -±-N(w-2)-^-N(w-l)+^-N(w)  =  0 

2AZ  AZ  2AZ 


X(w)  A(w)  B(w)  G(w) 

Between  the  boundary  conditions,  we  have  for  l<iz<w: 


'  N(iz  + 1)  -  2N{iz)  +  N{iz  - l)~ 

~N{iz  +  i)-N(iz-i)~ 

AZ2 

2AZ 

-A 


Rearranging,  we  have: 

N(iz+ 1)-^-  N(iz)+-^-  N{iz- 1)-—  N(iz+ 1)+-|-  A 
AZ2  *  AZ2  AZ2  2A Z  v  '  2AZ 

Grouping  like  terms  results  in: 


JVCi(fe)A(/z)=0 


;/z-l)-AJVCi(iz)A(iz)=0 
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3) 


1 

2AZ 


N(iz- 1) 


+ 


2d 
A Z2 


-^(t) 


A'(iz)  + 


AZ1 


1 

2AZ 


JV(iz +  1)  =  0 


A(iz)  B(iz)  D(iz) 

We  then  have  the  following  equations: 


G(iz) 


For  iz  = 

iz  =  1  B(l)  N(l)  +  D(l)  N(2)  +  X(l)  N(3)  =  G(l) 

iz  =  2  A(2)  N(l)  +  B(2)  N(2)  +  D(2)  N(3)  =  G(2) 

iz  =  3  A(3)  N(2)  +  B(3)  N(3)  +  D(3)  N(4)  =  G(3) 

iz  =  4  A(4)  N(3)  +  B(4)  N(4)  +  D(4)  N(5)  =  G(4) 

u  u 

u  u 

u  u 

iz=w-3  A(w-3)  N(w-4)  +  B(w-3)  N(w-3)  +  D(w-3)  N(w-2)  =  G(w-3) 

iz=w-2  A(w-2)  N(w-3)  +  B(w-2)  N(w-2)  +  D(w-2)  N(w-l)  =  G(w-2) 

iz=w-l  A(w-l)  N(w-2)  +  B(w-1)  N(w-1)  +  D(w-l)  N(w)  =  G(w-1) 

iz=w  X(w)  N(w-2)  +  A(w)  N(w-l)  +  B(w)  N(w)  =  G(w) 

We  can  then  represent  these  equations  by  a  matrix: 


B(l) 

D(l) 

X(l) 

0 

0 

=> 

A(2) 

B(2) 

D(2) 

0 

0 

0 

A(3) 

B(3) 

D(3) 

0 

=> 

0 

0 

A(4) 

B(4) 

D(4) 

=^> 

U 

u 

u 

u 

u 

u 

U 

u 

u 

D(w-3) 

0 

0 

B(w-2) 

D(w-2) 

0 

A(w-1) 

B(w-1) 

D(w-l) 

X(w) 

A(w) 

B(w) 

The  tri-diagonal  matrix  is  represented  as  follows: 
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m 

N(  2) 

N(  3) 

u 

A(w-2) 

A(w-l) 

iV(w) 


'5(1)  A(l)  +  +  X(Y)N(l) 

~G(  1)  ^ 

A(2)N(2)  +  B(2)N(2)  +  D(2)N(2) 

G(  2) 

A(3)N(3)  +  B(3)N(3)  +  D(3)N(3) 

G(  3) 

U 

II 

II 

II 

A(w  -  2 )N(w  -  2)  +  B(w- 2)N(w  -  2)  +  D(w  -  2N(w  -  2) 

G(w-  2) 

A(w  - 1 )N(w  - 1)  +  B(w  - 1  )N(w  - 1)  +  D(w  -  l)JV(w  - 1) 

G(w- 1) 

X(w)N(w)  +  A{w)  N(w)  +  B(w)N(w) 

G(w) 

The  Thomas  Method  may  now  be  employed  to  solve  the  tri-diagonal  matrix: 


N(l)  = 


m 

i- 


^(l)g(2) 

^(1)Z>(2) 

*m<2) 

B(l)D(2) 


m  mmv 

l  mm  c 


ETA(l) 


E(l) 


or 

N(l)  =  ETA(l)  -  E(1)N(2) 


Substitution  into  the  second  matrix  equation  results  in: 


N(2)  = 


G(2)-,4(2)£L4(1) 

B(2)-A(2)E(l) 


m 

B{2)-A(2)E(\) 


N(3) 


ETA(2) 


E(2) 


or 

N(2)  =  ETA(2)  -  E(2)N(3) 

The  following  can  be  done  for  every  iz  such  that 


- - N(iz+ 1 ) 

B(iz)  -  A(iz)E(iz  - 1) 


=  G(iz)-A(iz)ETA(iz-l) 

1  ;  B(iz)-A(iz)E(iz- 1) 

ETA(iz)  E(iz) 

Resulting  in 

N(w-2)  =  ETA(w-2)  -  E(w-2)N(w-l) 

N(w-l)  =  ETA(w-l)  -  E(w-l)N(w) 

The  last  expression  will  remain  the  same: 

X(w)N(w-2)  +  A(w)N(w-l)  +  B(w)N(w)  =  G(w) 

We  can  solve  the  last  3  equations  by  substitution,  resulting  in: 

_  G(w)  +  X(w)E(w  -  2)ETA(w- 1)  -  X(w)ETA(w  -  2)  -  A(w)ETA(w- 1) 
(W)  B(w )  -  A(w)E(w  - 1)  +  X(w)E(w  - 1  )E(w  -  2) 


We  can  now  solve  the  last  equation  as  it  only  has  one  unknown.  With  a  value  of  N  at 
iz  =  w,  we  can  now  solve  for  N  at  iz  =  w-1 .  Similarly,  a  value  can  be  solved  for  the 
previous  N  until  iz  =  1.  A  spreadsheet  such  as  Microsoft  Excel  is  most  useful  for  this 
exercise  and  is  enclosed  as  an  attachment.  After  a  value  for  N  at  each  iz  has  been 
obtained,  we  can  plot  N  versus  iz.  Additionally,  plotting  the  log  of  N  versus  iz  will 
indicate  an  overall  removal  efficiency. 
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